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Abstract 

We analyze the finite temperature chiral restoration transition of the (D = d + 1)- 
dimensional Gross-Neveu model for the case of a large number of flavors and fixed total 
fermion number. This leads to the study of the model with a nonzero imaginary chemical 
potential. In this formulation of the theory, we have obtained that, in the transition region, 
the model is described by a chiral conformal field theory where the concepts of dimensional 
reduction and universality do apply due to a transmutation of statistics which makes fermions 
act as if they were bosons, having zero energy. This result should be generic for theories with 
dynamical symmetry breaking, such as Quantum Chromodynamics. 
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1 Introduction 

Recently, Kocic and Kogut U have shown that, in the limit of a large number of flavors (N), 
the universality class of the finite temperature chiral symmetry restoration transition in the 3D 
Gross- Neveu model is mean field theory, in contrast to the "standard" sigma model scenario 
which predicts the 2D Ising model universality class. The responsible for the breakdown of the 
"standard" scenario, dimensional reduction plus universality Q, would be the absence of canonical 
scalar fields in the model, since the a auxiliary field has a composite nature. Indeed, as explained in 
both the density and the size of the a meson increase with the temperature in such a way that 
close to the restoration temperature the system is densely populated with overlapping composites. 
In opposition to what happens in the dynamics of BCS superconductors, here the constituent 
fermions are essential degrees of freedom, even in the scaling region. Thus, we conclude that, 
for this model, the effect of compactifying the Euclidean time direction is simply to regulate the 
infrared behavior and supress fluctuations. This fermionic model was first analyzed in || and it was 
further argued that the results would not be an artifact of the large N limit [ffl. (In addition, it was 
shown how the Ising point is recovered in 4D Yukawa models beyond the leading order in the 1/N 
expansion and why this should not happen in Gross- Neveu models.) Lattice simulations of this 
3D model have verified the predictions of the large N expansion at zero temperature, at nonzero 
temperature, and at nonzero chemical potential separatedly j5|. The results for the critical indices 
have been checked and improved by larger scale simulations enhanced by histogram methods 0. 
The conclusion was that the data are in excellent agreement with mean field theory and rule out 
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the Ising model values. However, the question concerning the nature of the transition for large but 
finite N in Gross-Neveu models remained open until the recent work of Reisz |?J where he suggests 
that the mean field exponents are indeed an artifact of the large N limit. By a combination of 
different techniques, like 1/N expansion, dimensional reduction and high-temperature series, it 
is shown how the 2D Ising model exponents are recovered beyond the leading order in the 1/N 
expansion in such a way that the "standard" scenario is still valid, if N is not large enough. 

In this letter we shall study further the large N limit of the Gross-Neveu model. We will show 
that, if we fix the particle number, it is possible to get non-trivial exponents even if N is infinitely 
large. This third set of exponents is consistent with the "standard" scenario with the difference 
that, as we will see, the reduced theory is not the linear sigma model but, instead, again a four 
fermi theory for massless fermions. According to the "standard" linear sigma model scenario, if 
dimensional reduction plus universality arguments hold, the finite temperature transition of the 
.D-dimensional Gross-Neveu model would lie in the same universality class of the (d = D — 1)- 
dimensional linear sigma model. This is not what is found and the reason why is easy to understand 
in the finite temperature Matsubara formalism: fermions do not have a zero frequency Matsubara 
mode. (The introduction of a real chemical potential would not change this results since it would 
simply cause a shift of the frequencies in the imaginary axis.) On the other hand, we will show 
that the imposition of a constraint which fixes the total fermion number gives rise to a fermionic 
zero frequency Matsubara mode. In this case, the transition has neither mean field nor scalar 
critical indices but rather chiral indices, a feature of theories with long range interactions. The 
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basic idea is to introduce a pure imaginary chemical potential in the partition function of the 
model. The imaginary chemical potential arises when we impose a constraint which fixes the 
particle number ||. This is generally done by inserting in the functional integral a functional delta 
function which inforces exactly the constraint. As a consequence, an auxiliary field is introduced 
in the problem. Usually, saddle-point solutions for this auxiliary field are imaginary implying in 
this way a real chemical potential. (This happens for example in the non-linear sigma model.) 
However, as we will see, real solutions are also allowed in the Gross-Neveu model with constrained 
particle number. Such a solution can be thought as an imaginary chemical potential. In this new 
"physical" system, we do have a zero frequency Matsubara mode associated to fermions affecting 
the infrared sensitiveness of thermodynamic quantities and, consequently, changing its critical 
indices. We observed that the new critical exponents are exactly the same as those corresponding 
to a zero temperature dimensionally reduced effective theory. In other words, after the decoupling 
of nonzero frequency Matsubara modes (dimensional reduction) the effective theory for the light 
degrees of freedom has a set of critical indices identical to that of the original zero temperature 
one (universality) , where one simply has to replace D = d + 1 by d. 
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2 The Gross-Neveu model at fixed total fermion number 

We start by considering a (D = d + l)-dimensional fermionic model with fixed total fermion 
number B described by the canonical partition function 

Z = Tr{e- /3H 5(N-B)}, (1) 

where f3 is the inverse of the temperature, H is some Hermitian Hamiltonian and 

N = J d^x^ty) (2) 

is the fermion number operator. We can rewrite the partition function as a functional integral 
over Grassmann fields 

Z = jf D^DtP 5 (y tT^ty) - B^j e~ A ^^ 

= [ D^D*{jD6e- A ^^ 9 \ (3) 
Jb.c. 

where we have used the integral representation of the functional delta function. In the above 
equation A is the action of the unconstrained model while A is given by 

r rt 3 - iB rt 3 

A = A+ d d x dr$i-a<$ + — drB. (4) 
J Jo p p Jo 

where the Euclidean gamma matrices are antihermitian and satisfy {7^,71,} = —25^. As usual, 
the Grassmann fields are integrated with antiperiodic boundary conditions (B.C.) with respect to 
the Matsubara time. Note that the Bose field 6 is a function of r but not of x. This is expected 



since 9 is an auxiliary field related to a global constraint with respect to space rather than a local 
one. 

We will use as fermionic action, the action for the massive Gross-Neveu model with N 
components. By rewriting it in its standard auxiliary field version we get the following Euclidean 
action for the fixed fermion number model 



S E (if>,il>,(T,0) = \l dr \ ^ x 



(d -> 9 \ 

-7°tt- + V + m + go + i— 7 ) -0 
dr (3 J 



(5) 



For bare fermion mass m — 0, there is a discrete chiral symmetry ip — ► 7 S 0>, ip —>■ —ipjs which is 
spontaneously broken whenever a non-vanishing condensate < tptp > is generated. The condensate 

(ip(x)ip(x)) oc J V9VaVi)Vi) ]fi{x)ij){x^ e - SB ^^°\ (6) 

serves as an order parameter of the transition. It is well known that this is the role played by 
the auxiliary field o in eq. fl5|). Thus, we will study the phase structure of the chiral symmetry 
restoration transition with the aid of the order parameter a. Note that eq. (|5[) is invariant by a 
partially local gauge transformation. By this we mean a gauge transformation which is local in 
the Matsubara time and global in the space. In fact, the transformation ip — > e^ip, ip — * e~ l< ^ip, 
leaves the action invariant provided that 9-^9 + /3d(fi/dr. Thus, 9 plays the role of a "gauge 
field". It should be observed that, due to the periodicity of the 9 with respect to the Matsubara 
time, the last term in eq. (JD is also gauge invariant: Jg drdcp/dr = <j>(0) — 0(0) = 0. Later we 
shall discuss the role of this gauge symmetry with respect to the zero mode of the fermion field. 
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3 Saddle-point evaluation of the functional integral 

Integrating out exactly the fermions in eq. @ we obtain the effective action 

(d -> 9 \ 

-7° 7^: + zf-V + m + ga + 

+ f fi dr( fd d xV + ^). (7) 



/o \J 2 /3 

Now we perform a saddle-point evaluation of the functional integral. From the structure of S e ff 
we know that the saddle-point solution is exact at large N. If N is large enough, the fluctuations 
of the u and 9 auxiliary fields average out and for this reason we look for uniform saddle-point 
solutions. By imposing stationarity with respect to a we obtain the gap equation 

V ^ f A d d q g 

where u n = ^j{n-\- ~), and we have used the fact that, to the leading order, the fermion self-energy 
S comes from the a auxiliary field tadpole: Y, = m — g 2 < x/ji/j >. Here A is an ultraviolet cutoff 
and we have defined p = 9/(3. To find the temperature dependence of the order parameter a one 
has to solve the gap equation near criticality. The critical temperature is determined by 

V ^ M d d q l 



II C/ 



with u; nc = 7p(n + |). On the other hand, stationarity with respect to 9 gives us 



^Tr[7 -SV(T-T';x-x')] = 5, (10) 
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where Sf denotes the fermion propagator given by the inverse of the argument of the log in Eq.(7). 
Eq.fllOl) is equivalent to the statement < ip'ij) >= B which fixes the mean number of fermions. 
Note that some care is needed to evaluate the trace in ([K]) since 7 • Sp is divergent for r = r' 
(for x = x' there is no problem due to the presence of the ultraviolet cutoff A). The same kind 
of difficulty is encountered in the theory of fermion systems, for instance Fermi liquid theory, 
when defining expectation values of particle number operators using the one-particle propagator 
r0|. In order to remedy this problem it is sufficient to evaluate 7 • Sp{r — r'; 0) with a temporal 



point-splitting r — r' = rj where 77 is a small positive number. After we let rj — > 0+ getting a 
finite result. This means that if we write the Fourier representation of 7 • Sf(t — t'; 0), the limit 
t — t' — > does not commute with the Matsubara sum. This is the standard procedure used 
in many-particle systems to evaluate expectation values of one-body operators fp~0f1 . Evaluating 
explicitly the trace in ([X0|) , we obtain 

1 ~ M d d q e iu>nV( p _ Qn ) 

d n hj (27r)dq2 + (p-a; n )2 + E2 4' 1 ] 

where 77 — > 0+ and we have introduced the barion density b = B/(NQ), Q being the volume. Note 
that the sum is divergent without the convergence factor e luJnV even if p = 0. It is divergent simply 



because it is not defined (like the series 1 — 1 + 1 — 1 + ...) [11 



The Matsubara sums in eqs.@ and ([O]) are done by contour integration and we obtain 

, f A d d a 1 / 1 1 \ 

y J 2tt ) d uo \e e e^ + 1 e^e-^ + ir V ; 
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for eq. (||) and 

^ d d q ( 1 



(13) 



(27r) d \e ie e^ + 1 e i6 e -P" + \ t 
for eq. (|TT| ) and we have defined to = i/q 2 + £ 2 . One saddle-point solution for # is 6> = ij3fi, fi 
playing the role of an ordinary chemical potential. This solution is non-singular in the infrared 
since it does not generate any zero mode. In this situation the exponents are mean field. However, 
the solution 9 = (2m + where m is an integer, is also possible if B is even. (We need B 
even in order to ensure that the free energy is a real quantity.) Note that every solution of this 
form corresponds to the same free energy and therefore we can choose any of them, say 9 = n. 
This particular solution is related to any other solution 9 = (2m + 1)tt by a gauge transformation. 
The existence of such solutions is a consequence of the partially local gauge symmetry we have 
discussed in the beginning of this paper. 



4 Scaling properties at fixed fermion number 

Let us discuss the implications of the solution 9 = n to the critical behavior of the model. First, 
by using eq. @ together with eq. (|8|) we can rewrite the gap equation for the order parameter a 



as 



TO V ~ /-A d d q (q 2 + (Pc - u nc f) (|) - (q 2 + (p - u n ) 2 + S 2 ) 

S + Pc n h J (2*y (q 2 + ( Pc - u nc f) (q 2 + (p - u n f + £ 2 ) " ' ' ' 

This form of the gap equation is particularly well suited for extracting critical indices since the 



problem reduces to the power counting of the infrared divergences in the integral over q [12|. For 
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6 = 7r we see that the above integral is infrared divergent as S — > 0. The critical indices are 
defined by (ipipj ~ t@, (ipip} t ~ m 1 / 5 , (S)| m _ +0 ~ t" etc and, since S ~< >? = ^ 
to the leading order. Above four dimensions the integral is infrared finite and the scaling is mean 
field. Below four dimensions, however, the S — > limit is singular and the integral scales as 
(^nJ- 2 + S 2 ) -5- . At a critical point t = away from the chiral limit we have m ~ Thus, we 

can easily see that, below four dimensions the f3 and S exponents are 

/? = — L_ * = d-l. (15) 

The remaining exponents are easily obtained as well, rj = 4 — d and 7 = 1, and one can check that 
they obey hyperscaling. If we compare these indices with those obtained from the zero temperature 
case: (3 = nz§> 5 = D — 1, r/ = 4 — D,j = 1 (for a recent review see [fL3l| ), we conclude that both 
sets define a chiral conformal field theory in d and (D = d+1) dimensions respectively. We should 
also note that the exponents flT5D are the same as for a zero temperature dimensionally reduced 
Gross- Neveu model in which the chiral symmetry restoration transition will occur as we approach 
a thermally renormalized coupling constant obtained during the reduction procedure. 

It should be noted that the real solution = n makes the model to behave as a bosonic model. 
Thus, this solution transmutate the statistics of the fermions into bosons. Unlike to statistical 
transmutation in Chern-Simons models, the present situation is not confined to d — 3. 
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5 Conclusions 



We have shown that, even in the absence of canonical scalar fields, we do obtain non mean field 
critical exponents for the finite temperature chiral restoration transition in a four-fermion theory 
provided we introduce a global constraint fixing the particle number. At the saddle-point level, 
it amounts to introducing an imaginary chemical potential which plays the role of a "gauge" 
field. This procedure gives rise to a zero frequency Matsubara mode which is interpreted, after 
the limit of high temperatures, as the only relevant degree of freedom of the reduced theory. 
The reduced theory is still a four-fermi theory at zero temperature and in a lower dimension 
whose parameters carry a temperature dependence resulting from the reduction procedure ITJJ . 
It also exhibits dynamical breaking of discrete chiral symmetry which is restored near the critical 
thermally renormalized coupling constant, in such a way that the universality class is the same as 
the original finite temperature model. 

This behavior should also be expected in any model with dynamical symmetry breaking, such 
as QCD, and at finite density, since, on general grounds, the (bulk) behavior of the system at real 
(fixed) chemical potential /i is identical to that at a fixed fermion number B, provided B is the 
mean fermion number for the system at chemical potential /i. 
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